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We aim to understand the role of Coulomb interactions as well as of different equations of state 
on the disappearance of transverse flow for various asymmetric reactions leading to same total mass. 
For the present study, the total mass of the system is kept constant {Atot ~ 152) and asymmetry 
of the reaction is varied between 0.2 and 0.7. We find that the contribution of mean-field at low 
incident energies is more for nearly symmetric systems, while the trend is opposite at higher incident 
energies. The Coulomb interactions as well as different equations of state are found to affect the 
balance energy significantly for larger asymmetric reactions. 

PACS numbers: 25.70.Pq, 25.70.-z, 24.10.Lx, 25.70.Mn, 21.65.Cd 



I. INTRODUCTION 

The heavy-ion physics has attracted much attention 
during the last three decades [H-Q- The behavior 
of nuclear matter under the extreme conditions of 
temperature, density, angular momentum etc., is a 
very important aspect of heavy-ion physics. One of the 
important quantity which has been used extensively to 
study this hot and dense nuclear matter is the collective 
transverse in-plane flow [J 0, @]. This quantity has a 
beauty of vanishing at a certain incident energy. This 
energy is dubbed as balance energy {Ebai) or the energy 
of vanishing flow (EVF) d, 0] . This beauty is due to the 
counterbalancing of attractive mean-field at low incident 
energies and repulsive nucleon-nucleon (NN) collisions 
at higher incident energies. The balance energy of 



the masses ranging from C + C to U 



t238 



U 



238 



at 



different colliding geometries was studied experimentally 
and theoretically and found to be sensitive with the 
composite mass of the system [a, |Ml as well as with the 
impact parameter of a reaction [a. lol. [Toj . 
With the passage of time, isospin degree of freedom 
in terms of symmetry energy and NN cross section is 
found to affect the balance energy or energy of vanishing 
flow and related phenomenon in heavy-ion collisions 

Ba[ii!,P!. 

Experimentally, Pak et al., studied the isospin effects 
on the collective flow and balance energy at central 
and peripheral collision geometries [llj. On the other 
hand, theoretically, this effect is studied by using the 
isospin-dependent Boltzmann Uehling-Uhlenbeck model 

(iBuu) n a [isj . and isospin-dependent quantum 

molecular dynamics (IQMD) model [1, [H, [H, I3 . 

As noted, balance energy is due to the counterbalancing 
of the attractive mean-field and repulsive nucleon- 
nucleon collisions. The Coulomb interaction in in- 



termediate energy heavy-ion collisions is expected to 
play a dominant role in balance energy due to its 
repulsive nature. These effects are supposed to be more 
pronounced in the presence of isospin effects [16]. The 
comparative study which will show the shift in balance 
energy due to Coulomb interactions in the presence 
of isospin effects by taking into account asymmetry of 
reaction in a controlled fashion is still missing in the 
literature. The second point is the asymmetry of the 
reaction. In some of the studies, the asymmetry of a 
reaction is taken into care, but not in other, which is 
very important to study the isospin effects [H, [T7| . 
The asymmetry of the reaction can be defined by the 
parameter 77 = | [At — Ap)/{At + Ap) |; where At 
and Ap are the masses of target and projectile. The 77 
= corresponds to the symmetric reactions, whereas, 
non-zero value of rj define different asymmetry of the 
reaction. It is worth mentioning that the reaction 
dynamics in a symmetric reaction (77 = 0) can be quit e 
different compared to asymmetric reaction [r] ^ 0) [l9|. 
This is due to the deposition of excitation energy in 
the form of compressional energy and thermal energy 
in symmetric and asymmetric reactions, respectively. 
The effect of the asymmetry of a reaction on the mul- 
tifragmentation is studied many times in the literature 
[igI IitI [19| . Unfortunately, very little study is available 
for the asymmetry of the reaction in terms of transverse 
in-plane flow. 

In this paper, we will perform the first ever study for the 
balance energy in terms of asymmetry of the reaction 
and then observe the effect of Coulomb interactions, 
symmetry energy, equations of state as well as different 
frame of references. The IQMD model used for the 
present analysis is explained in the Sec. -II. The results 
are presented in Sec. -Ill, leading to the conclusions in 
Sec.-IV. 
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II. THE MODEL 

The isospin-dependent quantum molecular dynamics 
(IQMD)[i^ model treats different charge states of 
nucleons, deltas and pions explicitly \2l\, as inherited 
from the Vlasov-Uehling-Uhlenbeck (VUU) model [2^ . 
The IQMD model was used successfully in analyzing 
the large number of observables from low to relativistic 
energies [H, [H, HH, [H. One of its version (QMD), 
has been very successful in explaining the subthreshold 
particle production (23| . multi-fragmentation [23, [25| . 
collective flow [1, l26l| . disappearance of flow [y], and 
density temperature reached in a reaction [IJ]- We shall 
not take relativistic effects into account, since in the 
energy domain we are interested, there is no relativistic 
effect [131 • The isospin degree of freedom enters into 
the calculations via both cross sections, mean field 
and Coulomb interactions The details about the 

elastic and inelastic cross sections for proton-proton and 
neutron-neutron collisions can be found in Refs. 15, 27]. 
In this model, baryons are represented by Gaussian- 
shaped density distributions 



1 -{T-r,(t)f -(p-piit:y)^.2L 

Mr,p,t)^—^e 2i e ^ . (1) 

Nucleons are initialized in a sphere with radius R — 
1.12^^/'^ fm, in accordance with the liquid drop model. 
Each nucleon occupies a volume of h!^ so that phase 
space is uniformly filled. The initial momenta are 
randomly chosen between and Fermi momentum pp. 
The nucleons of the target and projectile interact via two 
and three-body Skyrme forces and Yukawa potential. 
The isospin degrees of freedom is treated explicitly by 
employing a symmetry potential and explicit Coulomb 
forces between protons of the colliding target and pro- 
jectile. This helps in achieving the correct distribution 
of protons and neutrons within the nucleus. 
The hadrons propagate using Hamilton equations of 
motion: 



dfl d < H > dpi d< H > 

dt dpi ' dt dri 

with 

<i/>=<r> + <F>is the Hamiltonian. 

= E^ + EE / W,p.t)VHr'.r-) 

■i i j>i 

xfj{r',p',t)dfdr'dpdp'. (3) 

The baryon-baryon potential V^^ , in the above relation, 
reads as 



V ■'{r' r) — Vgf,y^^^ + Vy^ka^^a + Vcoul + ^Sym 

= hS{f' -r)+ t25[? - f)p'^-i(^^^) 

exp{\ r' -r\/ii) Z.ZjC^ 
(I r' — r I r' — r I 

+ u}-TlTi.b{r\~r-). (4) 

Po 

Where p. = OAfm, U = -&MMeV and ti = IQQMeV. 
The values of ti and <2 depends on the values of a, /3 
and 7 01 . Here Zi and Zj denote the charges of the 
i*'' and j"* baryon, and T^, Tg are their respective T3 
components (i.e. 1/2 for protons and -1/2 for neutrons). 

The parameters /i and ti, ,^4 are adjusted to the 

real part of the nucleonic optical potential. For the 
density dependence of the nucleon optical potential, 
standard Skyrme-type parameterizations is employed. 
The Skyrme energy density have been shown to be 
very successful at low incident energies where fusion is 
dominant channel [1^ [2^. The Yukawa term is quite 
similar to the surface energy coefficient used in the 
calculations of nuclear potential for fusion js^]. The 
choice of the equation of state (or compressibility) is 
still a controversial one. Many studies advocates softer 
matter, whereas, much more believe the matter to be 
harder in nature [l^] • We shall use both hard (H) and 
soft (S) equations of state that have compressibilities of 
380 and 200 MeV, respectively. 



III. RESULTS AND DISCUSSIONS 

As discussed earlier, asymmetry of a reaction is 
found to affect the phenomena of muti-fragmentation 
in intermediate energy heavy- ion collisions [13, [l9| . On 
the other hand, system mass dependence of balance 
energy is studied many times in the literature To 
check the effect of Coulomb interactions and asymmetry 
of a reaction on the balance energy, we have fixed 
{Atot = At + Ap — 152) and varied the asymmetry 
of the reaction just like this: 26Fe^^ +4,4 RvP^ (77 = 0.2), 

24Cr50 ^^102 ^ 3)^ ^^^^40 5^112 ^ 0.4), 

16532 -^50 5^120 (^ ^ 0.5), I45z28 +54 Xe^^^ (7? = 0.6), 

sO^^ -I-54 Xe^^^ {rj = 0.7). The asymmetry of a reaction 
with multi-fragmentation in this fashion is varied many 
times [13, [1^ • The whole reaction dynamics is studied 
at semi-central geometry by varying the incident energy 
between 50 and 250 MeV/nucleon with an increment 
of 50 MeV/nucleon by employing hard as well as soft 
equations of state. We have checked the stability of the 
reacting nuclei in laboratory (lab) as well as in center of 
mass (cm.) frame by taking into account the Coulomb 
interactions. Our main purpose here is to understand 
the effect of equations of state and Coulomb interactions 
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FIG. 1: Asymmetry dependence of directed flow in lab as 
weU as center of mass frame in L.H.S, while, R.H.S for the 
relative effect of Coulomb interactions. The different lines in 
the figure are representing the effect of symmetry energy and 
Coulomb interactions. 



on the energy of vanishing flow or alternatively, on the 
balance energy by taking into account the asymmetry of 
a reaction. 



The directed transverse flow is calculated using < P'^ 



where Y(i) and Px{i) are, respectively, the rapidity dis- 
tribution and transverse momentum of the z*'' particle. 

To check the effect of frame of reference, we display 
in Fig.l, variation of the asymmetry r/ on directed flow 
< P^*'' > in lab. as well as center of mass frame at 
incident energy of E = 50 MeV/nucleon. The top and 
bottom panels in the right hand side of the figure are 
representing the relative Coulomb effect with respect to 
the asymmetry of the reaction. This relative effect is 
calculated as: 



< AP; 



dir 



\Coul 



>~< P']^^ >Coul+Sym — < P^ 



. . (6) 

As is evident from the figure, directed fiow is found to 
increase in a systematic manner in CM. as well as in 



lab frame at E = 50 MeV/nucleon with asymmetry of 
the reaction. The inclusion of Coulomb interactions 
does not alter the conclusions. Note that the increase in 
the asymmetry is related with the increase in the N/Z 
ratio. Because the symmetry potential for the neutron 
rich systems is stronger compared to the neutron poor 
systems due to large relative neutron strength. Further- 
more, the symmetry potential is repulsive for neutrons 
and attractive for protons. On the other hand, more 
negative value of directed fiow (dominating the mean 
field) is observed in the absence of Coulomb interactions 
in center of mass as well as in lab. frames. This is due 
to the enhancement of the chemical and mechanical 
instability domains in the absence of Coulomb interac- 
tions . Similar type of study and conclusion was also 
performed for nuclear stopping in ref. [16]. 
Extensive study in the literature proves the stability 
of reactions in lab. frame, but, keep in mind that the 
reaction under consideration in these studies were sym- 
metric in nature. As is clear from the figure, asymmetric 
systems are found to be more stable in the center of 
mass frame compared to the lab frame. Moreover, if one 
consider lab frame, one is surely missing the effect of 
asymmetry. To further strengthen the stability of center 
of mass frame with asymmetry, the relative effect of 
Coulomb interactions < AP**" \coui> is studied in both 
frames. The relative effect of the Coulomb interactions 
is found to decrease with increase in the asymmetry 
of a reaction. The systematic decrease can be seen in 
center of mass frame with asymmetry, while very weak 
dependence of Coulomb interactions on the asymmetry 
is obtained in lab frame. For the further study, we 
have opted the center of mass frame because we want 
to see the shift in the balance energy due to Coulomb 
interactions, whose effect is clearly visible in center of 
mass frame. 

Before we proceed further, let us check the time evo- 
lution of directed transverse fiow. In Fig. [21 we display 
the time evolution of directed fiow from E = 50 to 200 
MeV/nucleon in center of mass frame for Soft (L.H.S) 
and Hard (R.H.S) equations of state. Note that the 
compressibilities of soft and hard equations of state are 
200 and 380 MeV, respectively 0. The time evolution 
is plotted in the absence of Coulomb interactions to see 
the maximum effect of asymmetry of a reaction on the 
directed in-plane fiow. The figure reveals the following 
points: 

a) . The quantity < P^**" > is observed to be constant 
throughout the whole distribution of time, while, the 
large variation is observed in the value at initial and 
final time steps when observed in the lab frame (6] 

b) . With the increase in the incident energy, the directed 
fiow is approaching towards more positive value. This 
is due to the well known fact of increase in the frequent 
NN collisions with increase in the incident energy [6|. 

>NoCoui+Syrrjr.'j^ r£]^Q bchavior of the directed flow with asymmetry 
of a reaction follows the opposite trend at E = 50 
MeV/nucleon as compared to other high incident ener- 
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FIG. 2: The time evolution of directed flow at difTerent in- 
cident energies in center of mass frame in tiie absence of 
Coulomb interactions. The left and right panels are for soft 
and hard equations of state, respectively. 



gies. It has been discussed many times in literature and 
also clear from the present findings that attractive mean- 
field is dominating at E = 50 MeV /nucleon coinpared 
to higher incident energies under consideration [ssj . 
This is due to the different mechanisms contributing at 
low and high incident energies within isospin-dependent 
quantum molecular dynamics. It was shown by us as 
well as by others [1, d, [l^l that symmetry potential 
dominates the physics at low incident energy, while, NN 
cross sections one major driven force at higher incident 
energies. Furthermore, at low incident energies, sym- 
metry potential is repulsive for neutrons and attractive 
for protons. With the increase in the asymmetry of a 
reaction, the number of neutrons increases and hence 
comparative repulsion due to neutrons also increases 
leading to less attractive value of the flow. It is also clear 
from ref. 17], that with the increase in asymmetry of 
a reaction, the participant zone decreases and spectator 
zone increases. At higher incident energies, where the 
effect of symmetry energy is negligible, the contribution 
of NN collisions comes from the participant zone whereas 
mean-field contribution comes from the spectator zone. 
Hence directed flow is found to be less positive with 
increase in the asymmetry of the reaction, 
d). The directed flow has less positive value with soft 



FIG. 3: Excitation function of directed flow at different 
asymmetries with and without Coulomb interactions at semi- 
central geometry. 



equation of state compared to hard equation of state. 
The less positive means the dominance of mean-field. 
This is due to the different compressibilities of hard (380 
MeV) and soft (200 MeV) equations of state. Naturally, 
more is the compressibility, more are the number of 
collisions and hence more positive is the directed fiow. 
This is indicating that the directed fiow is sensitive 
towards the equations of state. 

Finally, in Figs|3] and |31 we display the excitation 
function of directed fiow at different asymmetries from 
f] = 0.2 to 0.7. The value of abcissa at zero value of 
< P*'' > corresponds to the energy of vanishing fiow 
(EVF) or alternatively, the balance energy {Ehai)- The 
Fig. [3] shows the shift in the balance energy due to 
Coulomb interactions, while. Fig. |4] is representing the 
shift in the balance energy due to different equations 
of state. In Fig. |31 one sees a linear enhancement in 
the nuclear fiow with increase in the incident energy. 
This increase in the transverse fiow is sharp at smaller 
incident energies (upto 200 MeV/nucleon). If one goes 
to higher incident energies, the value gets saturated 
as discussed in ref [1]. We have displayed here the 
results upto 250 MeV/nucleon, since we are interested 
in and around balance energy. In the presence of 
Coulomb interactions, more positive value of the flow 
is obtained. This is due to the well known repulsive 
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FIG. 4: Excitation function of directed flow with hard and 
soft equations of state in the absence of Coulomb interactions 
for different asymmetries. 



nature of Coulomb interactions. At higher energies, 
the repulsion due to Coulomb interactions is stronger 
during the early phase of the reaction and transverse 
momentum increases sharply. The overall effect depends 
on the asymmetry of the reaction. If one looks at the 
balance energy, the shift in the incident energy towards 
the higher value is obtained at < P^"" >= with 
the asymmetry of the reaction. This is showing that 
with increase in asymmetry of the reaction and in the 
absence of Coulomb interactions, attractive mean-field 
is dominating the large region of incident energy. The 
systematics of the balance energy with asymmetry of 
the reaction is discussed in Fig. [S] 

As we have seen in Figl^l that different equations of 
state show sensitivity towards the directed flow with 
respect to the asymmetry of a reaction. The detailed 
analysis with soft (S) and hard (H) equations of state 
is displayed in Fig. The excitation function of 

directed flow follows similar trend as explained in Fig. 
[31 For nearly symmetric systems (77 = 0.2), the balance 
energy is found to be independent of the equations 
of state, however, resonable differences are observed 
at higher incident energies. More positive values of 
directed flow are obtained with hard equation of state 
compared to soft equation of state. This is true at 
all asymmetries from (77 = 0.3 to 0.7). This is due to 



FIG. 5: Power law dependence of balance energy with asym- 
metry of the reaction. The lower panel is representing the 
relative % effect of Coulomb interactions on the balance en- 
ergy. 



different compressibilities of hard (380 MeV) and soft 
(200 MeV) equations of state. With an increase in the 
asymmetry of a reaction, the shift in the balance energy 
towards the higher value of incident energy takes place 
with soft equation of state compared to hard one. This 
is consistent with the findings in the literature jTlj . 

To sum up, in FiglSJ we have displayed the asymmetry 
dependence of balance energy. We displayed the results 
for hard and soft equations of state by switching off the 
Coulomb interactions. In addition , for a comparative 
study, the results in the presence of Coulomb interactions 
with soft equation of state are also shown. All the lines 
are fitted with power law of the form Ebai — C{riy , 
where C and t are the constants. The values of r in 
the absence of Coulomb interactions for soft and hard 
equations of state are 0.375 and 0.282, respectively, 
while in the presence of Coulomb interactions for soft 
equation of state is r = 0.06067. 

If we compare the asymmetry dependence of balance 
energy with mass dependence, the trend is opposite 
0. It clearly indicates that if one wants to study the 
isospin effect, then one has to consider the systems 
with Atot = constant, otherwise one is not able to 
get the exact information about the isospin effects. It 
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is also clear from the figure, that shift in the balance 
energy is observed due to Coulomb interactions as well 
as due to equations of state with asymmetry of the 
reaction. The shift is more due to Coulomb interactions 
in comparison to equations of state, indicating the 
importance of Coulomb interactions in intermediate 
energy heavy-ion collisions. The higher balance energy 
is obtained with Coulomb-off + soft equation of state 
followed by Coulomb-ofF + hard equation of state and 
finally Coulomb-on + soft equation of state. 
For the further understanding, the relative percentage 
difference in the balance energy is plotted in the lower 
panel denoted by the quantity A.E^°i^^% given by 



IV. CONCLUSION 

Our present aim was to understand the influence of 
Coulomb interactions as well as equations of state on 
the dynamics of large asymmetric reactions in semi- 
central heavy-ion collisions. At low incident energies, 
the contribution of mean-field is more for the nearly 
symmetric systems, while at higher incident energies, 
opposite scenario is observed. The balance energy is 
found to increase with the increase in the asymmetry 
of the reaction. The balance energy is afi^ected by the 
Coulomb interactions compared to difli'erent equations of 
state. 
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The Ai?^j"'% is found to increase with the increase 
in the asymmetry of a reaction. This indicates shift of 
the nuclear matter towards the attractive mean-field 
region in the absence of Coulomb interactions with the 
asymmetry of the reaction. This difference A£'^°''' = 30 
MeV/nucleon at r] = 0.2, while it is 115 MeV/nucleon at 
r] = 0.7. The diff'erence of 90 MeV/nucleon in the shift 
of balance energy with asymmetry cannot be ignored. 
This is the first ever study and experiments are called to 
verify the results. 
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